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Abstract 

The Tropical Easterly Jet (TEJ) is a prominent atmospheric circulation feature 
observed during the Asian Summer Monsoon (ASM). The simulation of TEJ by the 
Community Atmosphere Model, version 3.1 (CAM-3.1) has been discussed in detail. 
Although the simulated TEJ replicates many observed features of the jet, the jet 
maximum is located too far to the west when compared to observation. Orography 
has minimal impact on the simulated TEJ hence indicating that latent heating is the 
crucial parameter. A series of aqua-planet experiments with increasing complexity 
was undertaken to understand the reasons for the extreme westward shift of the 
TEJ. 

The aqua-planet simulations show that a single heat source in the deep tropics 
is inadequate to explain the structure of the observed TEJ. Equatorial heating is 
necessary to impart a baroclinic structure and a realistic meridional structure. Jet 
zonal wind speeds are directly related to the magnitude of deep tropical heating. 
The location of peak zonal wind is influenced by off-equatorial heating which is 
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closest to it. Hence the presence of excess rainfall in Saudi Arabia has been shown 
to be the primary reason for the extreme westward shift of the TEJ maximum. 

Keywords: Tropical Easterly Jet Indian Summer Monsoon CAM-3.1 orography 
aqua-planet precipitation 


1 Introduction 


The Tropical Easterly Jet (TEJ) is one of the most defining features of the Indian Summer 
Monsoon (ISM) which itself is a part of the Asian summer monsoon (ASM). The jet 
is observed mostly during the ISM that is in the months of June to September. It 
has a maximum between 50°E-80°E, Equator-15°N and at ISOhPa. The TEJ was first 
discovered by Kpteswaram (119581) . It has a great influence on the rainfall in Africa 


(Hulme and Tosdevin 


19891) . The correct simulation of the TEJ is important for accurate 


seasonal predictions and weather forecasting. 

The Tropical Easterly Jet is believed to be influenced by the Tibetan Plateau. Previous 


studies (Flohn 


1968, 


Krishnamurtil Il9711) have highlighted the presence of a huge upper 


tropospheric anticyclone above the Tibetan Plateau in summer. The origin of the Ti¬ 
betan anticyclone itself has been attributed to the summertime insolation on the Tibetan 
Plateau. This sensible heating is widely reg arded as the primary reason for this region to 
act as an elevated heat source. lElohnI (11965^ first suggested that in summer, southern and 


southeastern Tibet, i.e. south of 34°N-35°N, act as an elevated heat source which changes 
the meridional temperature and pressure gradients and contribut es significantly to the 
reversal of high tropospheric flow during early June. iFlohnl (|l968l) showed that sensible 
heat source over the Tibetan Plateau as well as latent heat release due to monsoonal rains 
over central and eastern Himalayas and their southern approaches generates a warm core 
anticyclone in the upper troposphere at around 30°N. This was a primary mechanism 
for establishing the south Asian monsoonal circulation over south Asia. According to 
Koteswaram these winds are a part of the Tibetan Anticyclone which forms during the 
summer monsoon over South Asia. He believed that the southern flank of the anticy¬ 
clone preserved its angular momentum and b ecame an easterly at approximately 15°N. 
Contradicting Koteswaram, iRaghavanI (|1973^ opined that the upper tropospheric zonal 
component of the equatorward outflow from the Tibetan anticyclone does not agree with 
the law of conservation of angular momentum. However the thermal gradient balance 
was applicable and this was responsible for the origin of the TEJ. _ 


The TEJ is not simply a passive atmospheric phenomena. iHulme and Tosdevinl (119891) 
studied the impact of the TEJ on Sudan rainfall. El Nino events were suggested as a di¬ 
rect control over Sahelian rainfall via the TEJ. The decelerating limb of the TEJ showed 
interannual variations in location in relation to Eastern Pacific warming. Fluctuations 
in the jet were shown to be responsible for altering precipitation patterns in Sudan. 
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Camberlid (|l995fl also reported on the significant linkages between interannual varia¬ 


tions of summer rainfall in Ethiopia-Sudan reg ion and strength and latitudinal extent 


of upper-tropospheric easterlies. More recently 


Nicholson et al 


(I 2 OO 7 II showed that the 


wave activity of the TEJ influences the African easterly jet (A EJ). The AEJ in turn 
influences the weather patterns over the Atlantic coast of Africa. 


Bordoni and Schneider 


^|2008^ discussed the role of TEJ in modifying the monsoonal circulation. Upper-level 
easterlies shield the lower-level cross-equatorial monsoonal flow (which later becomes the 
Somali jet) from extratropical eddies. This made the angular momentum conservation 
principle applicable to overturning Hadley cell dynamics. In the larger picture this im¬ 
plies a strengthening of the monsoonal regime. This shows that the TEJ actively shapes 
the climate and weather pattern in regions not under the direct influence of the Indian 
Summer Monsoon. 

In fact the TEJ also shows intraseasonal and interdecadal variations in its location. 


Sathivamoorthv et al.l (120071) reported that the axis of the jet undergoes meridional move¬ 


ment in response to active and break phases of the Indian Monsoon. ISathivamoorthv 
^|2005^ also observed a reduction in the spatial extent of the TEJ between 1960-1990. Ac¬ 
cording to him the TEJ almost disappeared over the Atlantic and African regions. This 
reduction coincided with the 4-5 decade prolonged drought conditions over the Sahel re¬ 
gion. He speculated that these two phenomena were associated with each other. This 
explanation also implies that latent heating significantly influences the characteristics of 
the TEJ. 


Ya (|l98ll ) did laboratory experiments to simulate the heating effect of elevated land. 


He introduced heating in an ellipsoidal block resulting in vertical circulation, an anticy¬ 
clone in the upper layer and a cyclonic flow in the lower layer. The pattern w as found to 
be qua litatively similar to summer time atmospheric circulation in south Asia. iRaghavan 
(197^ discussed the importance of Tibetan Plateau and explained that the jet owed 
its existence the th e temp e ratur e gradients that was partly influenced by the Tibetan 
high. According to IWand (j2006l) the elevated heat source of the Tibetan Plateau was 
instrumental in prov i ding a n an chor to lo c ate th e Tibetan high. On the other hand 

(120071 1 have argued that orography plays a 


Hoskins and Rodwelll ( 1995 1 and 


Liu et al 


secon dary role in determining the position of the summertime upper tropospheric anticy¬ 
clone. Llingxi and Yihnil (jl989[ ) studied the TEJ at 200 hPa and fou nd that precipitation 


Chen and van Loon 


chang es in the west coast of India led to changes in the jet structure. 

( 1987|) observed that at 200 hPa the jet was weaker during El Nino and Indian drought 
events. The TEJ was weaker in the drought years of 1979, 1983 and 1987 and stronger 
in the excess mo nsoon years 1985 an d 1988. 

According to lZhang et al.l (|2002l) . the South Asian High had a bimodal structure. Its 
two modes - one the Tibetan mode and the other the Iranian mode, both of which were 
fairly regular in their occurrence, were mostly influenced by heating effects. The former 
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owed its existence to diabatic heating of the Tibetan Plateau, while the latter occurred 
due to adiabatic heating in the free atmosphere and diabatic heating near the surface. 

Thus it can be seen that while there have been attempts in the past to explain some 
aspects of the TEJ and its effects, the issue about the importance of orography on the 
TEJ has not yet been settled. In the previous discussions variations in the jet were in 
the presence of the Himalayas and Tibetan Plateau. Although the importance of latent 
heating on the TEJ is quite clear now, it is necessary to confirm the importance of 
orography in deciding the location of the TEJ. It is necessary to study the jet in a GCM 
and understand the relative roles played by latent heating and orography on the TEJ. 
An Atmospheric General Circulation Model (AGCM), CAM-3.1, has been used to study 
the impact of orography and monsoonal heat sources on the TEJ. The jet location and 
its response to have been studied by: 

(i) Modifying the deep-convective scheme and thereby changing the monsoonal heating 
pattern. 

(ii) Removing orography. 

A series of aqua-planet experiments have been conducted determine the factors influencing 
the location, structure and strength of the TEJ in CAM-3.1. 


2 The Tropical Easterly Jet in reanalysis 


For wind data, NCEP ([Kalnaj^^^^IJ, |199® and ERA40 (Ug^alu^e^^ 


been used, and GPCP ( Adler et all 12003 ) and CMAP ( Xie and Arkinl . Il997 ) data for 


2005) data have 


precipitation. Unless otherwise mentioned the time period chosen spans years 1980 to 
2002. The maximum zonal winds occur at 150hPa. From Fig. [T]it is seen that the TEJ 
peaks in July to August. Hence, in present work the focus is on the TEJ during the 
month of July when it first reaches its maximum value. 



Figure 1: Time series of 150 hPa zonal wind averaged between 50°E-80°E, Equator-15°N 
showing TEJ peaking in July-August (NCEP: 1980-2002 avg) 
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The magnitude of peak zonal wind and its location are listed in Table [TJ Due to the 
good agreement between zonal winds of NCEP and ERA40, only shown the former is 
shown. In Fig. I2al the zonal wind and location of peak zonal wind, henceforth Umax, 
for NCEP at 150hPa, is shown. The ‘cross-diamond’ indicates the horizontal location of 
Umax- The location of peak zonal wind of each individual year from 1980-2002 was found 
and then the means were calculated. The standard deviation for the zonal winds for 
NCEP and ERA40 was found to be 2 m s“^ and 1.7 m s“^ respectively. The meridional 
section of the zonal wind is shown in Fig. [Sa] This vertical cross-section corresponds to 
the longitude where the zonal wind is maximum. The meridional structure shows the jet 
peak lying between equator and ~20°N. A vertical equator to pole tilt towards higher 
pressure levels can be observed. The mean height of maximum easterly zonal wind speeds 
are higher at higher latitudes. In Fig. |4^ the velocity vectors and geopotential high is 
shown. The peak is not over the Tibetan Plateau but to the west of it. 

Table 1: Magnitude ( m s“^) and location of peak zonal wind, and centroid of mean 
precipitation 


Case 

reanalysis and observation 

Zonal wind 

Precipitation 

Peak 

Lon Lat Press 

Lon 

Lat 

NCEP 

ERA40 

GPCP 

CMAP 

34.48 

33.61 

65.1°E 10.3°N 150 

68.7°E 10.4°N 150 

96.0°E 

95.8°E 

11.9°N 

10.0°N 

Case 


CAM-3.1 simulations 

Zonal wind 

Precipitation 

Peak 

Lon Lat Press 

Lon 

Lat 

Ctrl 

noGlOrog 

47.19 

43.45 

42.5°E 8.8‘’N 125 

43.5°E 7.6‘^N 125 

76.7°E 

77.0°E 

9.7‘^N 

7.4°N 

Gtrl: default orography 
noGlOrog: no orography 


2.1 Impact of latent heating on the location of the TEJ 

The TEJ is not a stationary entity. Variations in its position in different years can also 
be observed. This is most strikingly observed in July 1988 and 2002. The location of 
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Figure 2: July horizontal zonal wind profile at pressure level where maximum zonal wind 
is attained, ‘cross-diamond’ is location of peak zonal wind: NCEP, Ctrl and noGlOrog. 
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(c) Ctrl 


(d) Ctrl 


100 



(e) noGlOrog 


(f) noGlOrog 


Figure 3: July cross-section of zonal wind profile at the longitude (left panel) and latitude 
(right panel) where Umax is attained (see Table[T]for details): NCEP, Ctrl and noGlOrog. 
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(a) NCEP 



(b) Ctrl 


Figure 4: July velocity vectors and geopotential height (xlO^m) at 150 hPa. ‘Cross¬ 
diamond’ is location of peak geopotential height. Contour for NCEP and Ctrl is 14.36 
and 14.43 respectively which is 99.5% of individual peak. 



the TEJ in during these two years is radically different. This is seen in Fig. [5^ where 
the 30 m s“^ zonal wind contour and Umax locations show significant differences. In fact 
the location of Umax is shifted eastwards by ^20° in 2002. These differences are seen in 
ERA40 data as well. The vertical structure in these two years is sim i lar to Fig. [3al 


It is well known that in India 2002 was a drought year (jSikka . 


2003, 


Bhat . 2006 ) 


while 1988 was an excess monsoon year. The relationship between precipitation and the 
jet location in the month of July for these two years is now clarified. Here in order to 
be self-consistent,both precipitation and zonal winds from NCEP data have been used. 
This is because reanalysis data is self-consistent and that the dynamic response of the 
atmosphere is entirely on account of forcing from the same data. If the TEJ is influenced 
largely by latent heating due to Indian Summer Monsoon then the July 2002 shift should 
be due to the mean rainfall being eastward shifted. This is clearly seen in Fig. I5bl where 
precipitation differences between July 1988 and 2002 have been plotted. The differences 
are quite striking. In July 2002 there is a clear eastward shift in rainfall with the maximum 
being in the Pacific warm pool and relatively little in the Indian region. This is in contrast 
with July 1988 where Indian region received high amounts of precipitation. The jet in 
1988 peaks at 60°E-65°E while in 2002 the maximum is in the southern Indian peninsula, 
implying a ^20° westward shift. The same behaviour was also found for ERA40 data 
and hence is not shown. 


3 The Tropical Easterly Jet in CAM-3.1 


3.1 Description of CAM-3.1 and experiment details 

The AGCM that has been used for the present work is the Community Atmosphere 
Model, version 3.1 (CAM-3.1). The finite-volume dynamical core using the recommended 
2° x2.5° grid resolution has been used for all simulations. T he time step is 30 minu t es and 
26 levels in the vertical are used. Deep conv e ction is the IZhang and McFarland (119951) 
scheme while shallow convection is the iHadJ ( 1994 ) scheme. Stratiform processes em¬ 


ploys the lRasch and KristianssonI (|l998l) scheme updated bv IZhang et all (120031) . Clou d 
fraction is computed using a generalization of the sche me introduced by ISlingd (jl989l) . 
The shortwave radiation scheme employed is describ e d in BrieglebI (jl992[ ) . The longwave 
radiation scheme is from iRamanathan and Downevi (|l986r) . Land surface fluxes of mo¬ 


mentum, sensible heat, and latent heat are calculated from Monin-Obukhov similarity 
theory applied to the surface. Climatological mean SST was specified as the boundary 
condition. Sea surface temperatures are the blended products that comb ine the global 


Hadl ey Centre Sea Ice and Sea Surface Temperature (HadlSST) dataset (IRavner et al 
2003 ) for years up to 1981 and Reynolds et al. ( 2002 ) dataset after 1981. 


The model was run in its default configuration for a five year period. This simulation 
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(a) NCEP zonal wind, 1988 (red), 2002 (blue) 



(b) NCEP precipitation difference: 1988-2002 


Figure 5: |(a)| 30 m zonal wind contour (dashed) in 1988 and 2002 showing shift in 
TEJ in response to heating, cross-sections are at Umax pressure level (150hPa), ‘cross¬ 
diamond’ shows the location of maximum zonal wind; |(b)| precipitation difference ( mm 
day“^), year 2002 subtracted from 1988. All data from NCEP. 
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is referred to as the control (Ctrl) simulation. Additionally another simulation (referred 
to as noGlOrog) has been conducted to check the influence of orography on the TEJ. This 
has also been run for five years with same boundary conditions but with orography all 
over the globe removed. This latter simulation is used to investigate the direct influence 
of topography on the TEJ. All the simulation results presented in this paper are based 
on five year means. 


3.2 Impact of orography on the simulated TEJ 


The 5 year average of maximum zonal wind and its corresponding location have been 
computed. The horizontal (Eigs. [2btcj) and meridional (Figs. IScleil profiles of the TEJ for 
Ctrl and noGlOrog simulations are shown. As with reanalysis, the vertical cross-section 
is at the location of zonal wind maximum. The existence of a Tropical Easterly Jet can 
be observed. In both the simulations the first noticeable feature is a ^30° westward shift 
of the simulated jet. This shift in the default Ctrl was also documented by 


Hurrell et al 


(|2006^ where they analyzed the 200hPa JJA zonal wind fields. The location of the peak 
zonal wind is virtually the same for Ctrl and noGlOrog simulations while the jet is weaker 
in the absence of orography. The zonal cross-section also shows that the zero line at 
500 hPa is sandwiched between the peaks in the TEJ and low-level Somali jet. One 
major difference appears in the Ctrl and noGlOrog simulation - the maximum in the 
Somali jet is greater and restricted east of ^50°E in the former. This was explained by 
Chakrabortv et al.l (|2008ll . The absence of orography in the latter caused the low-level 
wind to spread out and thereby reduce the intensity without compromising on the total 
flow. This is also seen in NCEP data (Fig. I3bll where the Somali jet is also maximum to 
the east of ~50°E. The meridional cross-section is very similar to real-planet and shows 
the familiar vertical equator to pole tilt with the maximum lying between equator and 
~20°N. The depth is also maximum on the poleward side. 

The impact of orography in determining the location and spatial structure of the jet 
in CAM-3.1 is thus minimal. The location of the peak zonal wind is virtually the same 
for Ctrl and noGlOrog simulations while the jet is weaker in the absence of orography. 
Although CAM-3.1 does show reasonable fidelity in determining the spatial features of the 
TEJ, the discrepancies in location and magnitude of the maximum velocity of the TEJ 
need to be understood. With the insight gained from section lOI the spatial distribution 
of simulated and observed rainfall is now studied. 


4 Spatial distribution of heating in reanalysis and Ctrl 

The westward shift of the simulated TEJ indicates that sensible heating from the Tibetan 
Plateau may not play a major role in the existence and location of TEJ. During the 
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monsoon season the major source of heating in the tropics is latent heating and hence it 
is necessary to look at the role of latent heating. 
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(a) Precipitation: GPCP 



(c) Precipitation: Ctrl 



(b) Precipitation: CMAP 



(d) Precipitation difference 


Figure 6: (a)[ (b)[ (c) July precipitation ( mm day ^). |(d)| July precipitation difference 
( mm day-^) (Ctrl-(GPCP+CMAP)/2). 


Figs. [6atel show the precipitation in the month of July of GPCP, CMAP and Ctrl. 
The difference between Ctrl and mean of GPCP and CMAP is shown in Fig. I6dl The 
observational data have been averaged since the precipitation patterns are very similar. 
The contrast between reanalysis and model simulations is quite striking. Most noticeable 
discrepancies in model simulations are (i) significantly reduced precipitation in northern 
Bay of Bengal, East Asia, western Pacific warm pool (ii) a significant precipitation tongue 
just south of the equator between 50°E-100°E and (iii) Spurious precipitation in the Saudi 
Arabian region which is quite prominent and equal in magnitude to the precipitation 
peaks in central Arabian Sea and south-western Bay of Bengal. This implies a major 
re alignment in the loca l heating pattern. This unrealistic precipitation has been discussed 


by 


Hiirrell et al 


( 120061 ). 


Thus from Table [T] and Fig. |6]it appears that westward shift in the peak of the sim¬ 
ulated precipitation is responsible for the TEJ to be centered east Africa. The strongest 
impact seems to be from the significantly high anomalous precipitation in the Saudi Ara¬ 
bian regi on which could p l ay a p art i n the extreme wes tward shift of the TEJ. 


As in 


Kucharski et al 


(20091) and 


Davis et al 


(|2012ri , precipitation is used as a proxy 


for latent heating. The centroid of precipitation (Pc) has been computed for in a region 
that is spatially quite significant and covers the monsoon region. Different regions for 
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reanalysis and CAM-3.1 simulations have been chosen since the precipitation pattern is 
spatially different. For reanalysis the west Pacific warm pool is also considered, while it 
is ignored for model simulations. The region chosen for reanalysis is 60°E-130°E, 16°S- 
36°N and for the simulations it is 40°E-110°E, 16°S-36°N (region demarcated by red 
boxes in Figs. mg). From Table [T] it is seen that the mean precipitation in simulations 
and reanalysis are in reasonable agreement although Ctrl experiment overestimated the 
precipitation by ~15%. The excess precipitation could be a cause for stronger jet speeds. 
The centroid of the precipitation, henceforth P^, is computed using the following equation: 


Pj^i 


( 1 ) 


where, 

Xc and Uc are the zonal and meridional coordinates of the Pc, 

Pi is the precipitation at each grid point, 

Xi and yi are the zonal and meridional distances from a fixed coordinate system, in each 
case the grid point where peak precipitation occurs. 

The distances are measured from a coordinate system centered at a point where the 
precipitation is maximum in the chosen region. There is a clear ~20° westward shift 
in the precipitation centroid in Ctrl with respect to reanalysis. Precipitation pattern of 
noGlOrog is similar to Ctrl and hence even in this case the shift is ~20° westward. The 
latitudinal differences are not significant. Choosing slightly different averaging regions 
does not significantly distort this relationship. 

Consequent to the westward shift in the precipitation centroid and TEJ, the peak 
geopotential height at 150hPa in the simulations has also shifted westwards by ^25°. 
This is clearly seen in Fig. |4] where the geopotential contour corresponds to 99.5% of 
the peak, the location of which is indicated by the cross-diamond. The velocity vectors 
clearly show the anti-cyclonic flow around the geopotential high. 

Though this offers evidence of the primacy of latent heating in determining the lo¬ 
cation of the jet, the validity of this argument is further demonstrated by showing the 
precipitation patterns of the CtrlTdi 2 and noG10rogrdi2 simulations. The mean precipi¬ 
tation and centroid are shown in Table [TJ Since the precipitation pattern is now similar 
to real-planet (not shown), the same region as in real-planet has been used to compute 
the location of the precipitation centroid. Figs. I7afbl show the precipitation and TEJ in 
these two simulations. The location of the TEJ agrees well with real-planet. Although 
the mean precipitation is more in comparison to the previous simulations, the peak zonal 
wind speed is almost the same. The geopotential high (Eig. in the Ctrlrdi 2 simula¬ 
tion is now almost same as that in NCEP (Fig. ISl) . The spatial separation between the 
precipitation centroid and Umax is included in Fig. [8l The zonal separation is also in 
closer agreement to real-planet and this is may be attributed to the increased similarity 
in spatial pattern of precipitation in real-planet and these two simulations. The similarity 
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in the precipitation patterns and TEJ structures in both these simulations as well as the 
previously discussed Ctrl and noGlOrog simulations suggests that it is the location and 
strength of the heat source and not orography that controls the TEJ. The only difference 
in precipitation between the CtrlT(ii 2 and noG 10 rogT£;i 2 simulations is that the former 
has additional heating to the north of 20°N around 90°E longitude. The absence of this 
peak in the latter further suggests that the Tibetan Plateau is not the primary reason to 
set up the temperature gradients that lead to the formation of the TEJ. Sensible heating 
due to the presence of Tibetan Plateau may be less important when dominated by latent 
heat release due to convective processes. 

Thus it has been demonstrated that there is a close correspondence between the source 
of heating, spatial location of TEJ and the geopotential high. The negligible influence 
of orograph y and strong effe ct of heating on upper level wind patterns has also been 


discussed by 


Liuetal 


(1200711 . It is also worthwhile to note that the Somali jet in NCEP 


and Ctrl simulation is roughly in the same location and hence the lower and upper-level 
winds (Figs. I3bld|) are not merely reflections of each ot her, not jus t in magnitude but also 
in location, as simple though insightful models, e.g. ICilll (Il980l) suggest. However the 
spatial location of the low-level jet and upper-level TEJ is more closely correlated in the 
noGlOrog simulation and this is due to the absence of orography as will also be shown in 
the subsequent discussion on aqua-planet simulations. 

However this mean pattern does not present the full picture. Between 70°E-100°E 
the precipitation peak in real-planet is more northwards in comparison to GAM-3.1 while 
the opposite is true between 30°E-60°E. In fact real-planet hardly shows rainfall in the 
latter region. Hence it is not clear how the spatial structure of the precipitation influences 
the positioning of the TEJ. This aspect will need to be clarified before further analysis 
and this will be the focus of the next section where a major simplification is adopted by 
running GAM-3.1 in the aqua-planet (AP) mode. 


5 Aqua-planet simulations 

Since orography has hardly any impact on the TEJ in the real-planet simulations it is 
instructive to study the atmospheric response only due to heating. The multiplicity of 
heat sources both in reanalysis and Ctrl preclude any easy interpretation of the influence 
that each heat source has on the TEJ. A major simplification is conceivable if one removes 
orography, land, sea-ice totally and further remove seasonal cycles, and yet retain all the 
important physics that the AGCM offers. This also implies that one needs to prescribe 
SSTs as a boundary condition. All this is possible in the aqua-planet configuration of 
CAM-3.1. The solar insolation is perpetually fixed at 21®* March which is March Equinox. 

In order to understand the role played by heating in determining the structure and 
location, CAM-3.1 has been run in aqua-planet configuration. The basic state consists 
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(a) Ctrlr(ii 2 (125 hPa) 



(b) noG10rogTdi2 (125 hPa) 


—30.m/s 



(c) CtrlTdi 2 (150 hPa) 


Figure 7: (a) (b) July precipitation ( mm day“^, shaded) and horizontal zonal wind ( m 
s“^) profile at pressure level where maximum zonal wind is attained; ‘cross-diamond’ is 
location of peak zonal wind, [(^ velocity vectors and geopotential height (xlO^m) at 150 
hPa; ‘cross-diamond’ is location of peak geopotential height; contour value is 14.46 which 
is 99.5% of peak: Cir\Tdi 2 and noG10rogT(ii2- 
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Figure 8: Separation (Pc-Ui) between location of precipitation centroid and maximum 
zonal wind. ‘Cross-diamond’ denotes aqua-planet simulations, ‘plus’ denotes real-planet 
simulations and reanalysis. 


of a uniform background SST on which additional heat sources are imposed. The heat 
sources are indirectly specified by setting an SST perturbation on this uniform SST back¬ 
ground. Precipitation induced on account of this SST perturbation is representative of 
total atmospheric heating. The implicit assumption is that latent heating is the dominant 
effect in the region of precipitation. 

The rationale for imposing heat sources on a uniform SST backgroun d in contrast to a 


zonall y symmetric and meridionally varying SST profile is now explained. iB aiendran et al 


(^2013^ showed that equatorial easterlies will be simulated even if a zonally symmetric but 
meridionally varying SST profile, symmetric about the equator, is used. The existence 
of an equato rial jet also depends on t he presence of twin or sing le Inter-Tropical Conver¬ 
gence Zone (IBaiendran et al.l (|201,‘lh , iNeale and Hoskina (j2001ll ). In such a meridionally 
varying SST profile it is more difficult to determine the role played by weak heat sources 
in the formation of a jet in aqua-planet. 

The aqua-planet simulations with just uniform background SSTs of 20°C and 25°C 
had a weak equatorial easterly of about 10-15 m s“^ peaking at ^250hPa. In this simple 
configuration, it is implicitly assumed that SSTs beyond 60‘^N/°S do not signif i cantly 
influence equatorial dynamics. This was also suggested bv iHoskins and Rodwelll (Il995ll 
when they conducted a series of experiments to understand the Asian summer monsoon. 
Thus a uniform background SST is deemed to be the simplest basic state on which 
additional heating may be imposed to study the influence of heating on the TEJ. 

The reasons for the existence of these easterlies is as follows: these simulations had 
a band of weak precipitation near the equator. This would naturally imply that moist 
parcels arise aloft implying the presence of a weak Hadley Cell. This means that air is 
mixed latitudinally. Therefore there must be westerly motion at higher latitudes and 
easterly motion at lower latitudes owing to angular momentum conservation. The time- 
averaged torque on the whole atmosphere due to surface friction must be zero, which 
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requires that there be both easterly and westerly winds. Thus easterlies must prevail 
near the equator. Drag on surface eas t erlies also transfers angular momentum from the 
surface to the atmosphere (| Schneider! (|2006l) '). The issu e of t he ex istence of equatorial 
upper tropospheric easterlies has also been discussed bv iLeel (Il999ll . According to her, 
in the deep tropics the horizontal transient eddy momentum flux accelerates the zonal 
mean zonal wind. Transient eddies of intraannual and interannual timescales were im- 
portant determining factors. Hadley cell dynamics was also important. According to 


Bordoni and Schneider! (|2008l) the Hadley cell approaches the angular momentum conser¬ 


vation limit because of these upper-level easterlies (upper-level easterlies strengthen when 
additional off-equatorial heating is present). 

The details of all the aqua-planet simulations are in Tables [5] and [3] In Fig. [31 the 
different locations where SSTs are imposed are shown. The names of the aqua-planet 
simulations start with ‘AP’. The magnitude of peak zonal wind and mean precipitation 
are listed in Table jS] The important point to note is that the SSTs have been chosen 
such that mean of precipitation exceeding 5 mm day“^ compares well with reanalysis 
and real-planet simulations. Fig. [5| shows the zonal and meridional distance between 
Pc and Umax- Equation o has been used to compute the centroid. The zonal wind 
for each month has not been individually computed and then averaged; the zonal find 
fields have been added and then averaged for the six month period. Hence these locations 
correspond to grid point values of the model. This method is acceptable since the heat 
source is stationary and the response too will average out over a six month time scale. 
The coarse model resolution implies that minor fluctuations in location will hardly distort 
the main observations and inferences. The vertical sections at the longitude and latitude 
where Umax is attained. Some representative cases are discussed. 



30"E 60"E 90"E 120”E 150"E 180" 


Figure 9: Schematic showing locations of imposed SSTs for aqua-planet simulations. 
Refer Table [2] for further details. (Not drawn to scale). 
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Table 2: Location and nomenclature of SST profiles for aqua-planet simulations. See 
table |3] for full nomenclature. Figure [9] shows the schematic of the locations. 


Region 

Shape of SST profile 

Location of SST peak 


Circle (C) 

20 ° diameter 


S 

Oval (Oa) 

50° major axis 

16° minor axis 

90°E,20°N 


Oval (Ob) 

90° major axis 
16° minor axis 


E 

Rectangle 

10 ° slope 

62°E,4°S to 94°E,2°S 

Circle 

20 ° diameter 

60°E,4°N 

B 



50°E,20°N 

C 

Circle 

20 ° diameter 

70°E,14°N 

D 



85°E,10°N 




120°E,20°N 

P 

Oval 

As in (Ob) above 

150°E,10°N 

130°E,2°S 


5.1 Single heat source simulations 

For single heat source simulations, the SST profiles imposed have circular and oval shapes 
all centered at 90°E,20°N which mimic the off-equatorial monsoonal heat source in the 
northern Bay of Bengal. The circular profiles (AP_S1 and AP_S2) have a diameter of 
20° while the oval shaped profiles have major axes of either 50° (AP_S3 and AP_S4) or 
90° (AP_S5 and AP_S6). All oval profiles have a minor axis of 16°. The simulations 
with oval-shaped SSTs serve to demonstrate the effect of the shape of the heating region 
on the jet. The peak SSTs are at 90°E,20°N. They go linearly down to the background 
temperature. However the profile with major axes of 90° (AP_S5 and AP_S6) has a non¬ 
linear SST gradient which is used to study if the gradient has any qualitative change on 
the jet. The background temperature is 25°C. In each set, the first simulation in each set 
has 29°C peak SST while the second has 32°C peak SST. Eor example, AP_S3 has 29°C 
and AP_S4 has 32° C SST peak on a uniform background temperature of 25° C. 

The meridional and zonal cross-sections (Figs. llOafbll are shown for AP_S6 simulation. 
Meridionally, the jet is not symmetric. The major difference with reanalysis and real- 
planet simulations lies in the vertical shape of the jet. With single off-equatorial heat 
source the jet has a pole to equator tilt towards higher pressure levels. This is opposite 
to reanalysis and real-planet (Figs. [3ate| . Below 500 hPa, there is also a region of weak 
low-level jet which is the aqua-planet counterpart of the Somali jet (not shown). Further, 
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Table 3: Magnitude ( m s“^) and location of peak zonal wind, and magnitude ( mm 
day“^) and centroid of mean precipitation for aqua-planet simulations. Refer Fig. [^and 
table [2] for shape and location of SST profiles. 


Single heat source (at region S) simulations 


SST Zonal wind Precipitation 



Peak 

Shape 

Peak 

Lon 

Lat 

Press 

Mean 

Lon 

Lat 

AP.S1 

29°C 

Circle (C) 

23.64 

70°E 

6 °N 

175 

8.81 

90.2°E 

21.9°N 

AP.S2 

32° C 

33.97 

80°E 

10°N 

150 

12.41 

90.5°E 

21.6°N 

AP.S3 

29°C 

Oval (Oa) 

30.75 

65°E 

6 °N 

175 

9.24 

88 .8°E 

20.7°N 

AP.S4 

32°C 

37.98 

67.5°E 

6 °N 

150 

14.74 

89.4°E 

21.4°N 

AP.S5 

29°C 

Oval (Ob) 

35.58 

45°E 

4°N 

175 

10.17 

86.3°E 

20.5°N 

AP.S6 

32°C 

42.56 

40°E 

4°N 

150 

15.06 

88 .6°E 

20.1°N 


Uniform background temperature: 25° C 


Multiple heat source simulations 


Case 

QQ'T 


Zonal wind 


Precipitation 


Peak 

Lon 

Lat 

Press 

Mean 

Lon 

Lat 

AP_E1 


15.54 

52.5°E 

2°N 

225 

11.52 

74.2°E 

3.6°S 

AP_E2 


17.58 

37.5°E 

6 °N 

225 

11.49 

72.5°E 

2.4°S 

AP.N1 

R 

33.08 

37.5°E 

10°N 

150 

9.02 

67.7°E 

16.7°N 

AP.N2 

0,0,JJ 

30.41 

57.5°E 

6 °N 

175 

9.27 

72.3°E 

16.2°N 

AP.NE1 

B,E 

26.73 

37.5°E 

8 °N 

175 

10.37 

67.0°E 

2.7°N 

AP.NE2 

23.80 

42.5°E 

14°N 

150 

7.51 

67.3°E 

3.1°N 

AP_M1 

R R 

30.89 

55°E 

6 °N 

125 

11.40 

70.8°E 

5.4°N 

AP_M2 


24.91 

52.5°E 

4°N 

150 

11.62 

71.6°E 

2.6°N 

AP_M3 


24.31 

40°E 

10°N 

150 

9.03 

71.9°E 

2.3°N 

AP_M4 

B,C,D,E 

36.02 

35°E 

8 °N 

150 

9.03 

68.7°E 

5.4°N 

AP_M5 


34.43 

32.5°E 

6 °N 

150 

10.86 

69.2°E 

4.9°N 

AP_M6 

B,C,D,E,P 

35.70 

32.5°E 

8 °N 

150 

10.22 

68 .8°E 

3.1°N 

Uniform background temperature: 22°C 


Peak temperatures: < 29°C 
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the low-level westerly was observed to be either very weak or almost non-existent for those 
simulations with reduced spatial extent of heating and/or lesser overall heating rates, for 
example AP_S2 and AP_S5. 

To the east of the jet, there is an upper-level westerly intrusion. A hint of this westerly 
is also observed in the reanalysis and real-planet simulations (Fig. El right panel). In Fig. 
EH the precipitation region is demarcated by the 5 mm day ^ and above shaded contours. 
It can be seen that precipitation (also refer TableEJ for the AP_S6 case has very high rates 
which is unrealistic. For comparatively low heating rates, a baroclinic structure similar 
to reanalysis and real-planet simulations does not exist. This only suggests that for such 
cases, the easterlies developed are more of an add-on to the equatorial easterlies developed 
with uniform background SSTs. However even low heating more than doubles the peak 
zonal winds generated in comparison to only uniform SST everywhere. This shows that 
even relatively higher zonal easterly wind speeds may not create vertical structures that 
have any resemblance with reanalysis and real-planet cases. So one cannot understand the 
structure of the TEJ with a single heat source. In all the single heat-source simulations 
the easterlies were present around the entire tropics. In each of the circular and oval sets, 
increased heating due to increased SSTs also increases the zonal wind speeds. Points A 
and B in Fig. El are for AP_S2 and AP_S6. 

The first thing to note is that compared to the change in mean precipitation mag¬ 
nitudes, the zonal separation differs significantly. The precipitation centroid also varies 
by less than ^5°. The maximum separation is for the 90° major axis oval-shaped SST 
(AP_S6) while the least is for the circular SST profile with 32°C peak (AP_S2). This 
demonstrates that more zonally constricted heating reduces the zonal separation between 
Pc and Umax- When the spatial extent of heating is less, the zonal wind has relatively 
rapid zonal acceleration and a slower zonal deceleration. Intense convergence in a rela¬ 
tively small region and subsequent south-easterly and then almost easterl y movement o f 
the air mass will cause zonally rapid acceleration and slow deceleration. iMishral (1198711 
attributed the southward movement of the jet to the spherical geometry of th e earth. 
Non-li nearity was also responsible for causing the jet to shift towards the equator (|Mishra 
(19931)). Increased spatial extent of heating causes increased acceleration length as is also 
evident from Table El It appears that spatially extensive heating imparts greater energy 
to the mean tropical flows, since higher zonal winds exist over longer zonal lengths. This 
is an important point and will be stressed upon in the discussion where multiple heat 
sources are introduced. The jet is thus to the south-west of the peak heating which is a 
departure from reanalysis and real-planet simulations. 

The location of the geopotential contour is very near to the region of maximum precip¬ 
itation. This is clearly seen in Fig. llOdI Between 70°E-II0°E the meridional velocities 
are also significant. This is a major departure from real-planet simulations where the 
contour is very much to the west of the region of precipitation. 
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(a) Meridional cross-section 



(b) Zonal cross-section 



(c) Horizontal cross-section: (150 hPa) 


30.m/s 



(d) Velocity vectors and geopotential height 


Figure 10: Zonal wind (ms profile. Vertical sections are at the |(a)| longitude and 
latitude where Umax is attained; |(c)| horizontal cross-section at pressure level where 
Umax is attained, ‘cross-diamond’ is location of Umax, ‘star’ is location of Pc, precipitation 
contours ( mm day“^) are shaded; [(dj] velocity vectors and geopotential height (xlO^m) 
at 150 hPa; ‘cross-diamond’ is location of peak geopotential height; contour value is 14.09 
which is 99.5% of peak: AP_S6 (Aqua-planet with single off-equatorial heat source in 
region S shown in Fig. 0. 
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Additional simulations were also done with SST peaks at 90°E,10°N. The uniform 
background temperatures were 20°C and 25°C. The magnitudes of the SST peaks ranged 
from 23°C to 37°C. For heating at 10°N, zonal wind speeds remained approximately 
~20 m s“^ and did not show any increase even with higher heating rates. In contrast, 
for heating at 20°N, zonal winds showed a monotonic increase with heating. The wind 
speeds increased from ~20 m s“^ to ^50 m s“^ as heating increased. Thus heating in 
lower latitudes does not result in increase in jet zonal velocities. 

5.2 Simulations with multiple heat sources resembling Ctrl sim¬ 
ulation 

The major finding of the previous section is that a single off-equatorial heat source is 
inadequate in explaining the structure of either the reanalysis TEJ or the full AGCM 
TEJ. Hence it is necessary to study the impact of combinations of heat sources. Fig. |9] 
shows the regions where precipitation exists in Ctrl simulation. An additional advantage 
that Ctrl simulation has more precipitating regions than reanalysis in the region of interest 
and so it is more beneficial to understand these heat sources in determining the location 
and strength of the TEJ. If these aqua-planet simulations have a jet structure resembling 
the TEJ in Ctrl or noClOrog simulation in location as well as in magnitude, then it is an 
important step forward. This is because it can be then claimed that the location of the 
TEJ as largely a construct of the location of heating and not on any complicated land- 
ocean interactions. This section is an attempt to understand the interactions between 
the various heating and how they interact to produce jet-like structures. 

The names AP_N1 to AP_M6 are for these cases. The locations and peak SSTs 
are given in Tables [21 [3] (also refer Fig. [9]). As explained before, in regions of overlap 
the maximum SSTs are chosen which removes any kinks in the profiles. The uniform 
background temperature is 22°C. Though it would have been more realistic to incorporate 
a 25°C background temperature, it was thought that 22°C would enable us to clearly 
decipher the role the different zones play as elsewhere precipitation would be less while 
at the same time 22°C background temperature would allow for minimal convection to 
occur. The peak SSTs again do not exceed 29°C and in most cases are less than that. 
In all cases, SSTs linearly increase to the peak value. Although only a few simulations 
of each type have been listed, more such simulations have been conducted with different 
SST peaks in the same locations in order to check the validity of the assertions made. 
They are not presented for want of space. 

5.2.1 Simulations without near-Equatorial heating 

These cases are AP_N1 and AP_N2 which have heating in regions B, C and D as shown in 
Fig. |9l This set has no equatorial heating. Since the heating is mostly off-equatorial, the 
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meridional and zonal structures are similar to those with single heat source at 90°E,20°N 
and hence are not shown. The increased heating caused a low-level westerly in both 
cases. The interesting point to note from Tableland Fig. [TT]is the ^25° westward shift 
of Umax in the 1®* simulation (Fig. Illal) in comparison to the 2"^^ simulation (Fig. Illal) . 
The zonal separation between and Umax (Fig. [5]) also shows a similar difference. The 
precipitation region in these two simulations shows that in the second case, precipitation 
region B is less than 10 mm day“^. This shows the importance of the heating that is 
closest to the peak zonal wind location. When the westernmost heating is below a certain 
threshold, the jet is influenced by the next closest heating, in this case the one in region 
C in Fig. m Simulations with heating only in regions C and D have also been conducted 
(not shown). In these cases, the zonal separation between Pc and Umax is about 15° 
less compared to AP_N1 and AP_N2. In other words, additional heating at location 
B increases this separation. This could be one of the reasons why Ctrl and noGlOrog 
simulations show a westward shift in the TEJ relative to reanalysis. In the real-planet 
simulations, there is significant precipitation in the Saudi Arabian region. 



(a) AP_N1 (150 hPa) 



(b) AP.N2 (175 hPa) 


Figure 11: Shift in location of Umax in response to locations of heating. Zonal wind ( 
m s“^) profile at pressure level where maximum zonal wind is attained. Precipitation 
contours ( mm day“^) are shaded, ‘cross-diamond’ indicates location of Umax, ‘star’ is 
location of Pc- Aqua-planet with multiple off-equatorial heat sources in regions B, C and 
D shown in Fig. [9l 


5.2.2 Importance of near-equatorial heating 

Referring to Fig. [6c] one can see that in the Ctrl simulation, during July, there is a broad 
precipitation zone just south of the Equator between 50°E-100°E and a zone just north of 
this zonal tongue centered at 60°E. Our goal is to induce similar precipitation patterns by 
imposing SST peaks in this region. Two such simulations have been conducted which are 
labelled as AP_E1 and AP_E2. The details are given in Tables |2l |3| and Fig. 0 Heating 
is present only in region E. Whenever there is an overlap, the maximum temperature is 
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used which ensures that there are no kinks in the SST profiles. However the background 
SST is now 22°C. The rationale for this choice will be explained in the next subsection 
where multiple-heat source simulations are discussed. Both rectangular and circular SST 
profiles have been imposed. The peak SSTs do not exceed 29°C and are made to increase 
linearly to the peak value. 




(a) AP_E2: Meridional section (b) AP_E2: Horizontal section: (225 hPa) 


Figure 12: Zonal wind (ms profile, (a) meridional section is at the longitude where 
Umax is attained; m horizontal cross-section at pressure level where Umax is attained, 
‘cross-diamond’ is location of Umax, ‘star’ is location of Pc, precipitation contours ( mm 
day“^) are shaded: Aqua-planet with multiple heat sources in region E shown in Fig. |9l 


From Table [3] it can be seen that the zonal wind speeds barely qualify as a jet. The 
speeds are hardly more than one-and-a-half times the easterly obtained for a uniform 
background simulation. The zonal wind structure is shown in Fig. 1121 The location of 
peak zonal wind is at a lower height (below 200 hPa). The most important observation is 
the stark contrast in the meridional structure vis-a-vis the off-equatorial heating. From 
the meridional section (Fig. ll2aL it is seen that the equator to pole tilt is towards higher 
pressure levels as in Ctrl and reanalysis. A low-level westerly of comparable strength 
extending below 500 hPa was also observed. Thus with just equatorial heating similar to 
real-planet, the vertical structure becomes baroclinic. Compared to reanalysis and real- 
planet there is a significant westerly to the east of the upper-level easterly. This signifies 
the reduced zonal extent of the easterly flows in comparison to single heat-source cases. 

But the horizontal structure (Fig. B2bll immediately shows that the zonal wind struc¬ 
ture is actually nowhere near the real-planet TEJ. The peak zonal wind was observed 
to be to the west of 90'”W and hence very far from Pc- This far-off value has not been 
tabulated since the region chosen for locating the peak zonal wind is 0-90‘’E. In spite of 
the maximum heating being in the south of the equator, the easterly is in the northern 
hemisphere. All these show that near-equatorial heating is necessary for imparting a baro¬ 
clinic structure somewhat resembling reality. But stand-alone near-equatorial heating is 
by itself insufficient to generate a jet structure that resembles the TEJ. 
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(a) AP_NE1: Meridional section (b) AP_NE1: Horizontal section: (175 liPa) 



(c) AP_NE2: Meridional section (d) AP_NE2: Horizontal section: (150 hPa) 


Figure 13: Zonal wind ( m s“^) profile. (a)|(c) meridional sections are at the longitude 
where Umax is attained; (b)|(d) horizontal cross-section at pressure level where Umax 
is attained, ‘cross-diamond’ is location of Umax, ‘star’ is location of Pc] precipitation 
contours ( mm day“^) are shaded: Aqua-planet with multiple heat sources in regions B 
and E shown in Fig. [HI 
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5.2.3 Interplay between heating in near-Eqnator and Saudi Arabian regions 

It is interesting to study the effect of heating in region B in addition to region E discussed 
above. These are the AP_NE1 and AP_NE2 simulations. It can be seen from Tables [1] 
and 131 that in these simulations the location of Umax is closer to Ctrl simulation. 

The meridional structure (Fig. I13all is similar to the equatorial heating case. There 
is also the equator to pole tilt that is observed in real-planet and real-planet cases. The 
meridional structures show the familiar equator to pole tilt (as in reanalysis) till the 
precipitation in region B (Saudi Arabia) does not increase beyond a threshold. This 
increase can be inferred from Fig. I13cl Even then the depth is moderated by the equatorial 
heating effects. The change in tilt is due to precipitation being more in region B as can 
be seen from Fig. I13dl in comparison to Fig. I13bl Here the equatorial heating in the 
latter is quite low and yet heating in B is not enough to increase the tilt in comparison 
to similar heating rates in the single heat source at 20°N). 

This shows that heating near the equator is essential in creating a meridional structure 
resembling reanalysis and real-planet. 

5.2.4 Effect of all heat sources except Pacific ocean warm pool region 

The discussion on the influence of and near-equatorial and off-equatorial heating can now 
be used to combine all of them and study how close the net effect is on the Ctrl simulation. 
These are the AP_M1 to AP_M5 simulations (refer Table [31). In this effort, once again 
the role of the easternmost heating on the jet structure and location is studied. The only 
difference between AP_M1 and AP_M2, and AP_M3 to AP_M5 simulations is the presence 
of heating in region B in the latter set. Heating in region D is expected to have minimal 
role to play as it is not the easternmost source. 

The first and foremost difference, as observed from Table |31 is in the zonal location of 
Umax- Precipitation in region B, as expected, causes a 15°-20° westward shift in the jet 
location while the location of Pc is hardly changed. This is clearly seen from Fig. |8]where 
points G and H are for AP_M1 and AP_M5 simulations. The zonal separation without 
heating in region B in the former causes the location of Umax to be much closer to Pc- 
The jet zonal velocities are generally greater and this is because of heating in regions B,C 
and D (see section [F.2.3|) . 

In all cases, the other zonal wind structures are similar in all cases (except for the 
difference in location of Umax when heating in region B is present). Hence in Fig. |T3| 
only AP_M5 simulation is shown. The upper-level meridional structure of the jet (Fig. 
[143 resembles real-planet simulations (Fig. |31 left panel). The equator to pole tilt is 
always towards higher pressure levels. This demonstrates that in aqua-planet configura¬ 
tion, precipitation patterns similar to real-planet simulations result in meridional shapes 
similar to full AGCM simulations. The zonal jet lengths are still less than real-planet 
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simulations.The eastward extent of the upper-level westerly intrusion is also unrealistic. 
This discrepancy will be addressed in section [S.2.51 below. 

The relatively weaker and more westward location of peak zonal wind speed in AP_M3 
simulation compared to AP_M5 is because the precipitation in region B in AP_M3 was 
much lower in comparison to AP_M5 (not shown). This once again underscores the 
influence of heating in region B (Saudi Arabia) in distorting the location of the jet in the 
Ctrl simulation. In spite of lower heating in region B in AP_M3 simulation, the meridional 
structure has the familiar equator to pole tilt (not shown). Except for the westerly to 
the east of 90°E, the overall TEJ simulated has a structure similar to Ctrl (Figs. [2b] and 
I3cldl) . 

The location of the geopotential high (Fig. I14dl) is now bears resemblance to Ctrl 
simulation (Fig. I4bl) . Although this peak in AP_M5 is shifted more eastwards in com¬ 
parison, from Tables |T] and |3| it is seen that the location of Umax as well as Pc is also 
westwards in comparison. On closer inspection, from Fig. |8| (points J and H) it can be 
observed that the zonal and meridional separation between Pc and Umax almost the same. 
Thus in the aqua-planet simulations if all the SST profiles had been shifted westwards by 
~10°, then the mosaic of precipitation, zonal wind and geopotential heights would have 
shifted westwards by about the same amount. Then there would be more similarity with 
real-planet simulations. 

It is also interesting to observe that the jet velocities increase by only ^5 m s“^ when 
compared with simulations with just an additional source at region B. In the previous 
single heat source simulations, the jet velocities increased when the intensity of heating 
increased both spatially and magnitude wise. Here, additional the heat source at 20°N 
does not significantly increase the jet strength. This shows that when multiple, well- 
distributed heat sources resembling reality are present, it is the intensity of heating and 
not the number of heating zones that determines the zonal wind strength. 

5.2.5 Importance of heating in the Pacific Ocean -warm pool 

In all the above multiple heat source simulations, the greatest mismatch in zonal wind 
structure with the real-planet simulations is the westerly that was always present to the 
west of 90°E. This westerly was prominently reflected in all the zonal cross-sections. 
From Fig. [6] it is observed that in reanalysis there is significant precipitation occurring 
east of 100°E. Although in this region, in both Ctrl and noGlOrog, the precipitation is 
less than observation, there was significant precipitation in excess of 5 mm day“^ in the 
Pacific warm pool which had not been incorporated in the aqua-planet configurations. 
Hence an additionale experiment, AP_M6, incorporating additional heating in region P 
as shown in Fig. |9l was conducted. The ratio of mean precipitation between 110°E-180° 
and 20°S-30°N for AP_M6 and Ctrl is ^0.95. Everywhere else the maximum SSTs were 
same as AP_M5. The results are shown in Fig. [121 The horizontal section (Fig. I15bp 
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(a) Meridional section 




30.m/s 



(d) Velocity vectors and geopotential height (150 
hPa) 


Figure 14: Zonal wind ( m s“^) profile. Vertical sections are at the | (a) | longitude and (a' 
latitude where Umax is attained; (c) horizontal cross-section at pressure level where Umax 
is attained, ‘cross-diamond’ is location of peak zonal wind, ‘star’ is location of pre¬ 
cipitation contours ( mm day“^) are shaded; [(d)] velocity vectors and geopotential height 
(xlO^m) at 150 hPa, ‘cross-diamond’ is location of peak geopotential height, contour 
value is 13.79 which is 99.5% of peak. AP_M5 (Aqua-planet with multiple heat sources 
in regions B, C, D and E shown in Fig. Ej). 
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now shows the extension of the easterly beyond 90°E and this is also reflected in the 
zonal section (Fig. IlSal) . Thus the acceleration length of the TEJ is now closer to real- 
planet and reanalysis. The meridional section is similar to AP_M5 simulation and hence 
is not shown. The location and magnitudes of peak zonal wind in AP_M5 and AP_M6 
are almost the same. Thus heating in the Pacific warm pool is necessary for the eastward 
extent of the TEJ. 




(b) Horizontal cross-section (150 hPa) 


Figure 15: Zonal wind ( m s profile, (a) zonal wind at the latitude where Umax is 
attained, m horizontal cross-section at pressure level where Umax is attained, ‘cross¬ 
diamond’ is location of Umax, ‘star’ is location of Pc, precipitation contours ( mm day“^) 
are shaded. AP_M6 (Aquaplanet with multiple heat sources in all regions, except S, 
shown in Fig. El). 


6 Conclusions 

The July climatological structure of the Tropical Easterly Jet in observations has been 
studied and compared with the simulations by an AGCM. The TEJ in reanalysis in July 
1998 and 2002 had significant zonal shifts. Reduced precipitation in the Indian region in 
2002 made the jet have its maxima in the southern Indian peninsula as compared to 1988 
when high rainfall in the Indian region resulted in the TEJ having its maxima over the 
Arabian sea region. 
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The TEJ in the Ctrl simulation has errors in the spatial location; otherwise the sim¬ 
ulated TEJ bears similarities with observations. Removing orography left the spatial 
location and structure is practically unchanged. This leads us to conclude that the TEJ 
is not directly influenced by orography. The primary reason for the shift in the simulated 
TEJ was because the location of precipitation in both Ctrl and noGlOrog is westwards 
when compared to reanalysis. 

Additional experiments were conducted to check if the TEJ is primarily influenced 
by latent heating. Changing the deep-convective relaxation time scale both in Ctrl and 
noGlOrog simulations confirmed this. In these new simulations the precipitation is more 
accurately simulated and most importantly anomalous precipitation in Saudi Arabia no 
longer occurred. The jet followed the shift in precipitation and relocated to the correct cli¬ 
matological position. The absence of orography once again had no impact on the location 
of the jet. This conclusively proves that the TEJ is independent of orography. Ghanging 
the default value of deep convective time-scale also demonstrated the secondary role of 
orography. Both Gtrlrdi 2 and noG10rogTdi2 have very similar precipitation patterns and 
hence the TEJ in both is correctly simulated. 

To understand why the TEJ was shifted westward in the AGGM, aqua-planet exper¬ 
iments were conducted. 

1 The simulation of TEJ in an aqua-planet configuration of the AGGM shows that orog¬ 
raphy and land-sea interactions are not as important as latent heat release. 

2 The total acceleration length is cirrelated to the zonal extent of the heating. 

3 A heat source at 20°N appears more to be robust in generating wind speeds that may 
be referred to as a jet while equatorial heating alone does not generate TEJ. 

4 Equatorial heating is necessary to generate a strong low-level westerly that imparts the 
vertical baroclinic structure to the TEJ. However it is insufficient in generating true 
TEJ horizontal structure. 

5 Equatorial heating is essential to create meridional structures seen in observations and 
full AGGM simulations. Greater poleward depth is possible only if equatorial heating 
is present. 

6 The longitudinal location of peak zonal wind is influenced by the off-equatorial heating 
that is closest to it. It has been demonstrated that rainfall in Saudi Arabia causes 
the extreme westward shift of the TEJ in full AGGM simulations in comparison to 
observations. 

7 Heating in the Pacific warm pool is essential to cause eastward extension and increased 
acceleration length of the TEJ. 
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When all the important heat sources are incorporated in the aqua-planet configuration, 
many observed features of the TEJ were simulated. Thus aqua-planet simulations play an 
important role in understanding the role of heat sources in the absence of any influence 
of land and orography. 
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